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We present an evolution equation for the Bjorken x dependence of diffractive dissociation on 
hadrons and nuclei at high energies. We extend the formulation of Kovchegov and Levin by relaxing 
the factorization assumption used there. The formulation is based on a technique used by Weigert to 
describe interjet energy flow. The method can be naturally extended to other exclusive observables. 
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QCD at very high parton densities is one of the most 
active frontiers both in high-energy and nuclear physics 
and one of the topics where both fields clearly profit from 
close collaboration. With the advent of the LHC in 2007 
this topic will further gain importance. Both the search 
for new physics in proton-proton collisions and the in- 
vestigation of high energy medium effects in heavy ion 
collisions require a solid understanding of multiple glu- 
onic interactions (at the very least for the analysis of 
backgrounds). Presently, the rapid output of precise ex- 
perimental data at RHIC, where the same effects should 
be present, though less pronounced, provides the main 
driving force behind new theoretical developments. One 
of the theoretically most attractive approaches is known 
under the name of color glass condensate (CGC) and 
one of its main elements is the JIMWLK-equation de- 
scribing the evolution of characteristic quantities with 
the squared cm energy s [2j. Our paper is based on this 
approach. 

At high energies, as reached in RHIC and LHC ex- 
periments, most QCD observables receive strong contri- 
butions from multiple soft gluon emission and multiple 
interactions of "hard" particles with soft gluons present 
in the event. A reliable and transparent method to re- 
sum the effects of these soft gluons on the hard lead- 
ing particles can be formulated by using gauge links 
U(x;y) — Pexp— ig J y dz.A(z) where the trajectories 
(from y to x) represent the quasiclassical paths of the 
hard particles while the soft gluons appear in the expo- 
nent. Previous work has focused on inclusive reactions. 
Here we demonstrate how to extend this program to ex- 
clusive reactions and work out the example of diffractive 
dissociation, where we can compare to a known limiting 
case that emerges if we use a factorization assump- 
tion as in the reduction of the JIMWLK to the Balitsky- 
Kovchegov (BK) |j,|5| equation. 

As an example let us recall that e.g. the total cross 
section of deep inelastic scattering of a virtual photon on 
a nuclear target can be written in terms of the Us as: 

vms(Y, Q 2 ) = I da [d 2 r \yj 2 \(r 2 a(l - a)Q 2 ) 

Jo J (1) 

ld 2 b (tr[l - U x Ul]/N c ) Y 



where \ip 2 \(r 2 a(l — a)Q 2 ) is the probability of a photon 
to split into a quark- antiquark pair of size r = x — y 1 
carrying longitudinal momentum fractions a and 1 — a, 
respectively. The remaining integral over the impact pa- 
rameter b = (x + y)/2 yields the cross section of a qq 
dipole of size r. The rapidity Y = In 1/x is taken to be 
large. The gauge links U x and t/t represent the leading 
hard quark and antiquark within the virtual photon wave 
function. They propagate at fixed transverse coordinates 
x and y along straight lines from = -oo to = oo: 



Pexp— ig I dz b + (z ,x,x + = 0) 



(2) 



In J5J we have anticipated that the hard particles inter- 
act, to leading order, only with 6 + , the soft component 
of the gauge field with rapidities below some Yq . Generi- 
cally A»(x) = U l (x) +5A> Jl (x) with b+ = 8{x )f3{x) and 
b~ = b = 0; 5A^ denotes all hard fluctuations. 

The averaging indicated in (JJJ is over the soft fields 
b + and represents all the interactions with the target 
through gluons softer than the original qq. As such it 
contains nonperturbative information that can not di- 
rectly be calculated. Since this decomposition into hard 
and soft modes is rapidity dependent, the dipole cross 
section turns rapidity dependent as well. By considering 
hard corrections 6 A one can systematically calculate the 
Y dependence and find RG equations for the dipole cross 
section. A direct approach leads to an infinite hierarchy 
of equations, the Balitsky hierarchy Q, which can only 
be solved after truncation. A more compact formulation 
in terms of a single diffusion equation can be given in a 
functional language. To this end one parametrizes the 
lack of knowledge about the averaging procedure by us- 
ing a functional Zy [U] , which takes on the meaning of a 
statistical distribution function: 



) Y ■= / D[U]...Z Y [U] 



(3) 



D[U] is a Haar measure that is normalized to 1. The s-, 
respectively Y-, evolution for the dipole cross section and 
all the other more complicated correlators in the Balitsky 
hierarchy is then given by the JIMWLK equation, which 
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governs the evolution of the functional weight Z: 



d Y Z[U] Y = -H. 



JIMWLK 



Z[Uh 



(4) 



describes a Fokker-Planck type diffusion problem in 
a functional context. The JIMWLK Hamiltonian 



H 



JIMWLK 



2tt 2 ' 



(•5) 



(integration over repeated transverse coordinates is im- 
plied here and below) contains a real emission part 
proportional to a new adjoint Wilson line U£ b that 
signals the appearance of a new gluon in the final 
state and virtual corrections that guarantee finiteness 
of the evolution equation. The remaining ingredients 
are the kernel JC xzy = (x-z) 2 lz-y)' 2 an< ^ functional 
derivatives iV° that respect the group valued nature 
of the [/-fields: := — \Uxt a ]ijd / 8U x ,ij corresponds 

to the left invariant vector field on the group mani- 
fold, while its right invariant counterpart is given by 
i^l := [PUjijS/SUvv = -U^iVl 

To prepare for the treatment of non-inclusive observ- 
ables, we will now sketch how to recover this evolution 
equation from the underlying real emission amplitudes. 
The treatment parallels that of jet observables in |6(. To 
facilitate this construction, let us introduce Wilson lines 
Uy, x that are (slightly) tilted w.r.t. the lightcone. The 
derivation is then based on the observation that the whole 
cloud of Y ordered real gluons accompanying any given 
number of hard partons that can be characterized as a 
product of Wilson lines U^ ] xi ■ ■ ■ U$ can be gener- 
ated by the application of a single operator 

U[U,£] = Py 2 exp [»/ ' dY 1 dY 2 9(Y 1 -Y 2 ) 



ji jjab pb,i -™ 
J xz u Y 2 ,z^Y 2 ,z l v Yi ,x 



(6) 



with J* := ~r~r t~~ — ~tt the eikonal current in trans- 
verse coordinate space. The £ fields represent the glu- 
onic final states. The derivatives now act as iV Yx := 
— [UY,xt a ]ijS / SU~Y,x,ij; F-ordering is such that the hard- 
est gluon is rightmost. This ensures that gluons can only 
emit softer ones. For DIS, the hard seed consists of a qq 
pair represented by a product of Wilson lines U Y ^ x U Yy 
at projectile rapidities. [The external color indices are 
those of the amplitude.] Diagrammatically we get 



follow the logarithmically enhanced contributions. The 
remaining gluons indicate soft interactions with the tar- 
get below Y , which build up the U fields. They corre- 
spond to the initial condition of the evolution process. To 
recover the real emission part of the dipole cross section, 
we need to square this amplitude and integrate over phase 
space for the resolved gluons. The average over the soft, 
unresolved gluons is done separately in amplitude and 
complex conjugate amplitude Q: we distinguish corre- 
sponding eikonal factors U and U. The average over the 
resolved final states can be made explicit by averaging 
over the final state variables £ with a Gaussian weight. 
For an arbitrary functional F[£] this is expressed as |6j 



(8) 



The ^-correlator M^^ v = appropriately 
normalized, is given by 

MgJgv = ^S ab ^5^) v 6 YuY2 e(Y-Y 1 )e(Y 1 Y ) . (9) 

M is diagonal in coordinates and rapidities and restricted 
to the resolved evolution interval [Yq , Y] . In the exponent 
of (JHJ integration and summation over all indices is un- 
derstood. The resolved contribution of real emissions to 
the dipole cross section is thus obtained by the Gaussian 
average (JSJ over the functional 



G real K] = U[£/,e]U[E/,£] 



tr[l - (UW)y.x(UU^I 



Y,y\ 



N r . 



(10) 



where N£ y = tr[l - {UW)y, x {UW) Y J/N c is the dipole 
operator of the total cross-section. No matter how the av- 
erage over the non-resolved modes below Yq is achieved, 
the evolution of the complete real emission part is deter- 
mined by the Y dependence of the resolved contributions: 



= (u[u,mu,^jc, 



— ± -2- M S— 

2 iil 7X1 
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U Y y y KyK ■ 

(11) 



In this equation everything outside the shower operators 
only contains U factors at the upper Y limit. This allows 
to recast both of these averages in terms of averages over 
Wilson lines at this highest rapidity: one can set 
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# of gluons & 
Jlowcd insertion 




(7) 



where the vertical dashed line denotes the final state, 
where each line ends in a factor £. These stand for explic- 
itly resolved partons in an interval [Yq, Y] over which we 



{...} Y - (U[U, Z]U[U, i, so ft= J D[U]D[U]...Z Y [U,U] 

(12) 

We may drop the now unnecessary F-label on the Wilson 
lines. This result, as all inclusive quantities, only depends 
on products Uffi in the hard operators appearing in l|ll|) 
and thus also in the weight Z. By a redefinition UW — > 
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U, the integration over U then reduces to a factor of 1 
and one is back at © . then leads to the contribution 
of real emissions to the evolution equation 



{dyZ[U}) 



real 



\K** V U? iV^iVt J[U] 



(13) 



Inserting virtual corrections by the requirement of real- 
virtual cancellation in absence of interaction, one recov- 
ers the JIMWLK evolution as stated in Q ■ Exclusive 
quantities on the other hand will depend separately on 
U and U and require to keep both fields in Z along with 
more complicated evolution equations. 

Since exclusive quantities are characterized by specific 
restrictions on the phase space of produced gluons the 
physically most transparent derivation of a corresponding 
evolution equation is built on a systematic construction 
of the contributing real emission amplitudes. The first 
modification clearly concerns the ^-correlator M used to 
implement the phase space integrals. Diffractive disso- 
ciation, which corresponds to a rapidity gap on the side 
of the target, requires a factor u(k) = 9{Yk — Vgap) for 
each final state gluon with momentum k . (The gap ra- 
pidity Yg a p is assumed to lie in the resolved range.) The 
major change, however, results from the appearance of 
additional diagrams that disappear in the inclusive re- 
sult through complete real virtual cancellation. While 
for JIMWLK it is sufficient to consider branching pro- 
cesses that occur before the interaction with the Lorentz 
contracted target, exclusive observables like diffractive 
dissociation will receive contributions from reabsorption 
and production in the final state, i.e. after the inter- 
action with the target as shown in Fig^ Reabsorption 



St, 
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FIG. 1: Generic diagrams for exclusive processes with final 
state interactions. In the diagrams, rapidity of gluons in- 
creases both vertically, in the final state, and horizontally, 
with the distance of their emission vertex to the target: To 
leading logarithmic accuracy, ordering in Y coincides with 
ordering in z~ towards the interaction region. Consequently, 
emissions into the final state after the interaction do not it- 
erate: lines marked in the graph to the right are suppressed. 

of a gluon after the interaction in the amplitude takes 
a form similar to a virtual correction in the JIMWLK 
case, but contains the soft interaction with the target, 
i.e. a factor U per hard particle. Technically, the neces- 
sary diagrams can be constructed by introducing a "three 
time formalism" in which we distinguish z~ = — oo, = 
and = +oo as the times at which the initial hard parti- 
cles are created, the interaction takes place and the final 
state is formed respectively. The transition amplitude 



from z~ = — oo to +oo is then created in two steps: we 
use a shower operator to create gluons before the interac- 
tion but anticipate that some of them directly reach the 
final state while others will be reabsorbed after the inter- 
action. In order to also generate the final state contribu- 
tions with a shower operator, we introduce an auxiliary 
Gaussian "noise" S with the same average and correla- 
tor as in © and JjJJ . Furthermore we artificially split the 
U factors of the interaction region into two Wilson lines 
W and according to U = FWL (One may think of 
them as Wilson lines extending over the intervals [— oo, 0] 
and [0, oo], respectively, they will disappear in the final 
result.) We then obtain the full set of diagrams: 



pjf[s.e]£ 




(14) 



£ 




where the sum is over the number of gluons and allowed 
insertions. The dashed line through the interaction re- 
gion represents the auxiliary split of the Wilson lines into 
W and with accompanying 5 factors. The shower op- 
erators are given by 



Ui[E,£] = Pr 2 exp i / dY 1 dY 2 e(Y 1 -Y 2 )4 



U f [H,£] = Py 2 ex P i dYxdY 2 e(Y x -Y 2 )J. 



(15a) 



(15b) 



Eventually combining the above expression for the ampli- 
tude with the corresponding expression for the complex 
conjugate amplitude and differentiating w.r.t. Y yields 
all real emission contributions to the evolution Hamil- 
tonian as well as the interacting virtual ones. One still 
misses virtual lines that do not cross the interaction re- 
gions. These are again reconstructed on the level of the 
evolution equation. We obtain the full Hamiltonian: 

H = u(k)H r +H v + Hy (16) 

where the real gluonic corrections are produced by 



The remaining terms correspond to virtual corrections in 
amplitude and complex conjugate amplitude respectively 



2tt 2 
-a s 
2^ 
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(The last terms in these expressions are the interact- 
ing parts.) The evolution equation parallels J3J, with 
Z replaced by Z[U, U]. Note that H v and H„ taken in- 
dividually have the form of the JIMWLK-Hamiltonian: 
H v is the evolution Hamiltonian for the dipole opera- 
tor of the forward amplitude N° y = tr[l - U X U£\/N C , 
which, via the optical theorem, determines the evolu- 
tion of the total cross-section. Real contributions only 
occur outside the gap, as mandated by the factor u(k). 
If we remove that restriction by setting u(k) = 1, we 
expect complete cancellation of final state contributions 
and again a reduction to JIMWLK. Indeed, setting u(k) 
to 1 and acting with (|To|) on the dipole operator = 
tr[l - (UW) x (UW)l]/N c (which depends only on prod- 
ucts (Uffl)W), we find that the evolution Hamiltonian 
(|16|l reduces to the JIMWLK-Hamiltonian for Wilson 
lines (UW). The average over N xy then is written in 
terms of Z[U, U] = Z\UU^] with evolution according to 



dyZ[uW] = -H[UU% 1W m^Z[UW] 



(17) 



cancellation is complete. 

The operators that replace tr[l — U x U y ]/N c (i.e. N F ) 
in Q for diffractive dissociation of a photon are different 
in- and outside the gapjthe structures closely resemble 
the factorized results of |3(. In the gap, no colored object 
enters the final state, thus the initial qq, if in the gap, 
interacts with the target and emerges in a singlet state: 
here N xy , the operator for cross sections with rapidity 
gaps larger than Y gapi is a product of two traces: 



N D = N° -N° = 

xy xy yx 



tr(l - U x Ul) tr(l - UlUy) 



Nr. 



(18) 



it takes the form of the "square" of two "elastic" dipole 
operators corresponding to the two amplitude factors. 
Here, evolution of amplitude and complex conjugate are 
completely uncorrelated ([H v ,Hy] = 0, H r is absent): 



Z Y [U,U] = 



_ -(H v +H e ){Y-Y ) 



Z Yo [U,U), (19) 



and one may factorize Zy[U,U) — > Zy[U]Zy[U] unless 
the initial condition contradicts this. 

For Y > Y gap , production of gluons in the final state is 
allowed and the gg-pair can appear also in an octet state. 
Adding the octet part to 118JI removes one trace: 



N£ v = tr[(l - UlU x ) (1 - UiU v )]/N c 

= n°, + n?„ - iv£, . 



(20) 



The second line exposes further structure: With the ini- 
tial conditions on evolution for (N D ) imposed by (|18fl . we 
find that (N ) acquires the interpretation of the cross 
section of events with rapidity gaps smaller than Y gap . 
Following our previous reasoning we conclude that the 
average over the three operators in the second line of 
(|2*U|| can be described by using functionals Z[U], Z[U] 



and Z[UW], with their respective evolution given by a 
JIMWLK-Hamiltonian. Even if additional structure in 
the initial conditions does not prevent these simplifica- 
tions, initial conditions for the individual terms are dif- 
ferent from each other (c.f. (|18|l % ) and the inclusive case. 

The relation to the results of Kovchegov and Levin Q 
parallels the reduction step from JIMWLK to BK: There 
one observes that JIMWLK evolution of S xy [U] = 1 — 
N xy = tr(U x lI y )/N c takes the simple form 



dY{S xy )Y 



2tt 2 



d zlC XZ y(S xz S Z y S X y)Y , (21) 



Factorizing {S xz S zy ) 



where K xzy = (iB _ z)a(z _ y)a ■ 

{S xz ){S Z y) truncates the infinite Balitsky hierarchy and 
leaves us with the BK equation. For Y > Y gav where 

N£y = 1 - S X y[U] - Sy X [U] + S X y[UU^} Silld B&Ch 

S xy [..] obeys H21jl. the same reasoning leads to the evo- 
lution equation presented as Eq. (9) in [2|. lfl%|l and 
(fH3|l imply the required initial condition {N xy )(Y gap ) — 

(N xy ) 2 {Y gap ) for evolution above the gap. 

To summarize: We have developed a method which 
allows to generalize the JIMWLK approach to a large 
class of exclusive observables, by simply adapting the 
phase space constraints. We have worked out the ex- 
ample of diffractive dissociation. For all generalizations 
it is crucial to start from IR safe observables, otherwise 
reconstruction of virtual contribution via real virtual can- 
cellations must fail. 
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